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CN ; Abstract 

^ ■ The current observational and experimental bounds on time variation 

^ ' of the constants of Nature are briefly reviewed. 

^ : 1 Introduction 

<^ : 1.1 History 
00 ' 

To our knowledge, Dirac appears to have been the first who argued for the possibility 
of time variation of the constants of nature. As is well-known, dimensionless numbers 
involving G are huge (or minuscule). For example, the ratio of the electrostatic force to 
the gravitational force between an electron and a proton is Q 

^. ^ 

where e is the electric charge, rup is the proton mass and mg is the electron mass. Similarly, 
the ratio of the Hubble horizon radius of the Universe, Hq^ to the classical radius of an 
^ ■ electron is ^ 

; N2 = - 3 X 10^°/i-\ (2) 

where h is the Hubble parameter in units of kms~^Mpc^^. Curiously, the two nearly 
coincides, which motivated Dirac to postulate the so-called the large number hypothesis 
|Jl[|. In his article entitled "A new basis for cosmology", he describes 

Any two of the very large dimensionless numbers occurring in Nature are con- 
nected by a simple mathematical relation, in which the coefficients are of the 
order of magnitude unity. 

Thus if the (almost) equality Ni = C(l) x N2 holds always, then G must decrease with 
time G oc or the fine structure constant, a = e^, must increase with time a oc t^/^ 

since H oct'\ 

Nowadays we know that such a huge dimensionless number like A^^i is related to the 
gauge hierarchy problem. In fact, the gauge couplings are running (however, only log- 
arithmically) as the energy grows, and all the gauge couplings are beheved to unify at 
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the fundamental energy scale (probably string scale). The fact that A^i nearly coincides 
with N2 may be just accidental, and pursuing the relation between them is numerological 
speculation (or requires anthropic arguments). However, Pandora's box was opened. 



1.2 Modern Motivation 

String theory is the most promising approach to unify all the fundamental forces in nature. 
It is believed that in string theory all the coupling constants and parameters (except 
the string tension) in nature are derived quantities and are determined by the vacuum 
expectation values of the dilaton and moduli. 

On the other hand, we know that the Universe is expanding. Then it is no wonder to 
imagine the possibility of time variation of the constants of nature during the evolution 
of the Universe. 

In fact, it is argued that the effective potentials of dilaton or moduli induced by 
nonperturbative effects may exhibit runaway structure; they asymptote zero for the weak 
coupling limit where dilaton becomes minus infinity or internal radius becomes infinity 
and symmetries are restored in the limit 0, Thus it is expected that as these fields 
vary, the natural "constants" may change in time and moreover the violation of the weak 
equivalence principle may be induced p, |[ (see also 0, ||] for earlier discussion). 

In this article, we review the current experimental (laboratory, astrophysical and geo- 
physical) constraints on time variation of the constants of nature. In particular, we 
consider a, G, and A. We assume Hq = 100/ikm/s/Mpc with h = 0.65 for the Hubble 
parameter and Qm = 0.3 and = 0.7 for the cosmological parameters. For earlier 
expositions, see & |10[ for example. 



z a 

In this section, we review the experimental constraints on time variation of the fine struc- 
ture constant. The results are summarized in Table 1. 



2.1 Oklo Natural Reactor and a 

In 1972, the French CEO (Commissariat a I'Energie Atomique) discovered ancient natural 
nuclear reactors in the ore body of the Oklo uranium mine in Gabon, West Africa. It 
is called the Oklo phenomenon. The reactor operated about 1.8 billion years ago corre- 
sponding to 2; ~ 0.13 for the assumed cosmology {h = 0.65, Qm = 0.3, = 0.7). 

Shlyakhter noticed the extremely low resonance energy {E^ = 97.3meV) of the reaction 

149^^ + ^^ 150^^ + ^, (3) 

and hence the abundance of ^^^Sm (one of the nuclear fission products of "^^^U) observed 
at the Oklo can be a good probe of the variability of the coupling constants [^. The 
isotope ratio of ^^^S'm/ ^^^5"^ is 0.02 rather than 0.9 as in natural samarium due to the 
neutron flux onto ^'^^Sm during the uranium fission. From an analysis of nuclear and 
geochemical data, the operating conditions of the reactor was inferred and the thermally 
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averaged neutron-absorption cross section could be estimated. The nuclear Coulomb 
energy is of order Vq ~ IMeV, and its change AVq is related to the change of a as 



Vq Vq a ' 

By estimating the uncertainty in the resonance energy, Shlyakhter obtained a/ a = 
lO^^V^'^^- Damour and Dyson reanalysed the data and obtained (—6.7 ~ 5.0) x lO^^^yr^^ 
12 1 . Using new samples that were carefully collected to minimize natural contamination 



and also on a careful temperature estimate of the reactors, Fujii et al. reached a tighter 
boundQ (-0.2 ± 0.8) x lO'^V"^ II- 



2.2 (Hyper)Fine splitting and a 

Since fine structure levels depend multiplicitely on a, wavelength spectra of cosmologi- 
cally distant quasars provide a natural laboratory for investigating the time variability 
of a. Narrow lines in quasar spectra are produced by absorption of radiation in inter- 
vening clouds of gas, many of which are enriched with heavy elements. For example, the 
separation between the wavelength corresponding to the transition "^81/2 — P3/2 and 
the wavelength corresponding to the transition ^5*1/2 —^^ P1/2 in an alkaline ion is pro- 
portional to a^. Because quasar spectra contain doublet absorption lines at a number of 
redshifts, it is possible to check for time variation in a simply by looking for changes in the 
doublet separation of alkaline-type ions with one outer electron as a function of redshift 



1^ , [151 . By looking at Si IV doublet, Cowie and Songaila obtained the limit up to 2; ~ 3: 
|Aq;/q;| < 3.5 X 10"'^ [|I6|. Also by comparing the hyperfine 21 cm HI transition with 
optical atomic transitions in the same cloud at z ~ 1.8, they obtained a bound on the 
fractional change in a up to redshift z ~ 1.8: Aa/a = (3.5 ± 5.5) x 10~^, corresponding 
to a/a = (—3.3 ± 5.2) x lO^^^yr^^ |16|. Recently, by comparing the absorption by the HI 



21 cm hyperfine transition (at z = 0.25, 0.68) with the absorption by molecular rotational 
transitions, Carilli et al. obtained a bound: |Aa/a| < 1.7 x 10~^ [0. 



Webb et al. ||T^ introduced a new technique (called many-multiplet method) that 
compares the absorption wavelengths of magnesium and iron atoms in the same absorbing 
cloud, which is far more sensitive than the alkaline-doublet method. They observed a 
number of intergalactic clouds at redshifts from 0.5 to 1.6. For the entire sample they 
find Aa/a = (— 1.1±0.4) x 10~^, consistent with a null result within 3 cr. They noted that 
the deviation is dominated by measurements at ^ > 1, where Aa/a = (— 1.9±0.5) x 10~^. 

Recently, Webb et al. [|19| pesented further evidence for time variation of a by re- 



analysing the previous data and including new data of Keck/HIRES absorption systems. 
The results indicate a smaller a in the past and the optical sample shows a 4 cr deviation 
for 0.5 < ^ < 3.5: Aa/a = (-0.72 ± 0.18) x 10"^ They noted that the potentially 
significant systemtic effects only make the deviation siginificant. If the result is correct, it 
would have profound implications for our understanding of fundamental physics. So the 
claim needs to be verified independently by other observations. 



''They noted that data is also consistent with a non-null result: (—4.9 ± 0.4) x lO^^^yr^^, indicating 
an apparent evidence for the time variabihty. However, from the analysis of the isotope compositions of 
Gd, the consistency of the Sm and Gd results supports the null results. 
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2.3 Laboratory Tests: Clock Comparison 



Laboratory constraints are based on clock comparisons with ultrastable oscillatiors of 
different physical makeup such as the superconducting cavity oscillator vs the cesium 
hyperfine clock transition [^] or the Mg fine structure transition vs the cesium hyper- 
fine clock transition |2l[]. Since a hyperfine splitting is a function of Za, such a clock 
comparison can be a probe of time variation of a. The clock comparisons place a limit 
on present day variation of a and are repeatable and hence are complementary to the 
geophysical or cosmological constraints. Comparisons of rates between clocks based on 
hyperfine transitions in alkali atoms with different atomic number Z (H-maser and Hg"*" 
clocks) over 140 days yield a bound: \ala\ < 3.7 x 10~^''yr~ 



2.4 Cosmology and a: Big-Bang Nucleosynthesis and Cosmic 
Microwave Background 

Big-Bang Nucleosynthesis. The '^He abundance is primarily determined by the 
neutron-to-proton ratio prior to nucleosynthesis which before the freeze-out of the weak 
interaction rates at a temperature Tf ~ IMeV, is given approximately by the equilibrium 
condition: 

(n/p)~exp(-g/T;), (5) 

where Q = 1.29MeV is the mass difference between neutron and proton. Assuming that 
all neutrons are incorporated into ^Hg, ^Hg abundance Yp is given by ~ 2{n/p)/[l + 
[n/p)]. The freeze-out temperature is determined by the competition between the weak 
interaction rates and the expansion rate of the Universe. The dependence of Tf on the 
weak and gravitational couplings is given by 

Tf oc (6) 

where Gp is the Fermi constant. 

Changes in Yp are induced by changes in Tf and Q. However, it is found that Yp is 



most sensitive to changes in Q p3[. The a dependence of Q can be written as [24, UB, 



Q ~ 1.29 - 0.76 X Aa/a MeV. (7) 

Comparing with the observed Yp, a bound on Aa/a is obtained: |Aa/a| < 2.6 x 10~^ 
p5| . Recently, it is argued that the presently unclear observational situation concerning 



the primordial abundances precludes a better bound than |Aa/a| < 2 x 10 ^ p7| . 



Cosmic Microwave Background. Changing a alters the ionization history of the 
universe and hence affects the spectrum of cosmic microwave background fluctuations: it 
changes the Thomson scattering cross section, cjr = 87ra^/3mg, and hence the differential 
optical depth f of photons due to Thomson scattering through f = XeUpar where Xe 
is the ionization fraction and Up is the number density of electrons; it also changes the 
recombination of hydrogen. 

The last scattering surface is defined by the peak of the visibility function, g(z) = 
e~'^^^''dt/dz, which measures the differential probability that a photon last scattered at 
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redshift 




d/a(yr i) 


Atomic Clock m 







< 3.7 X 10-14 


Oklo ( D amour- Dyson ) 


0.13 


(-0.9 ~ 1.2) X 10-^ 


(-6.7-5.0) X 10-1^ 


Oklo(Fujii et al.|13[) 


0.13 


(-0.36 ~ 1.44) X 10-s 


(-0.2 ±0.8) X 10-1^ 


HI 21 cm[|T|] 


1.8 


(3.5 ±5.5) X lO-'^ 


(-3.3 ±5.2) X 10-1^ 


HI 21 cmH 


0.25,0.68 


< 1.7 X 10-5 




QSO absorption line[18 


0.5 - 1.6 


(-1.1 ±0.4) X 10-5 




QSO absorption line 19 


0.5-3.5 


(-0.72 ±0.18) X 10-5 




CMB|2g 


10^ 


< 10-2 ~ 10-3 


< 10-12 ~ 10-13 


BBN|23] 


1010 


< 2 X 10-2 


< 1.4 X 10-12 



Table 1: Summary of the experimental bounds on time variation of the fine structure 
constant. Aa/a = (athcn — «now)/«now A bound on Aa/a for CMB is a possible bound. 



redshift z. As explained in increasing a affects the visibility function g\z): it increases 
the redshift of the last scattering surface and decreases the thickness of the last scattering 
surface. This is because the increase in a shifts gi^z) to higher redshift since the equilibrium 
ionization fraction, xf'*^, is shifted to higer redshift and because more closely tracks 
rj,EQ larger a. 

An increase in a changes the spectrum of CMB fluctuations: the peak positions in 
the spectrum shift to higher values of i (that is, a smaller angle) and the values of Cg 
increase p8[ . The former effect is due to the increase of the redshift of the last scattering 
surface, while the latter is due to a larger early integrated Sachs- Wolfe effect because of 
an earlier recombination. It is concluded that the future of cosmic microwave background 
experiment (MAP, Planck) could be sensitive to |AQ;/a| ~ 10-2 — 10-^ [H, |29||.p| 



3 G 

In this section, we review the experimental constraints on time variation of the Newton 
constant. For more detailed review see ||33|. The results are summarized in Table 2. 



3.1 Viking Radar-Ranging to Mars, Lunar-Laser-Ranging and 

G 

If we write the effective gravitational constant G as G = Go ± G'o(t — to), the effect of 
changing G is readily seen through the change in the equation of motion: 

d2x _ GMx _ GoMx Go G^M x(t - to) .g^ 

dt'^ r3 r3 Go r 

^After this talk, new analyses of BOOMERanG [|o) and MAXIMA and the first year results 
from DASI |3|] appeared. In those data the second (and even the third) peak is now clearly seen, 
which discourages the motivation for non-standard recombination scenarios that were attracted attention 
regarding the interpretation of the apparent absence of the second peak in the previous data. 
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redshift 


AG/G 


G/G(yr-i) 


Viking Lander Ranging [p4| 







(2 ±4) X 10^12 


Lunar Laser Ranging |36| 







(1±8) X 10-12 


Double Neutron Star Binary [p8| 







(1.10 ± 1.07) X 10-11 


Pulsar- White Dwarf Binary |39|| 







(-9± 18) X 10-12 


Helioseismology |43| 







< 1.6 X 10^12 


Neutron Star MassH 


- 3 ~ 4 




(-0.6 ±2.0) X 10-12 


BBN|6| 




-0.3 ~ 0.4 


(-2.7 ~ 2.1) X 10-11 



Table 2: Summary of the experimental bounds on time variation of the gravitational 

constant. AG/G = (Gthen - Gnow)/G'now 

Thus time variation of G induces an acceleration term in addition to the usual Newtonian 
and relativistic ones, which would affect the motion of bodies, such as planets and binary 
pulsar. 

A relative distance between the Earth and Mars was accurately measured by taking 
thousands of range measurements between tracking stations of the Deep Space Network 
and Viking landers on Mars. From a least-squares fit of the parameters of the solar 
system model to the data taken from various range measurements including those by 
Viking landers to Mars (from July 1976 to July 1982), a bound on G is obtained: G/G = 
(2 ±4) X 10-1 V"^ H- 

Similarly, Lunar-Laser- Ranging measurements have been used to accurately determine 
parameters of the solar system, in particular the Earth-Moon separation. From the anal- 
ysis of the data from 1969 to 1990, a bound is obtained: G/G = (0.1 ± 10.4) x lO^i^yr-i 



|3^; while from the data from 1970 to 1994, G/G = (1 ± 8) x lO^iV"^ 



3.2 Binary Pulsar and G 

The timing of the orbital dynamics of binary pulsars provides a new test of time variation 
of G. To the Newtonian order, the orbital period of a two-body system is given by 

'~ [Gm) ~ G2m2(l - e^f/^ ' ^ ' 

where a is the semi-major axis, i = r'^cf) is the angular momentum per unit mass, m 
is a Newtonian-order mass parameter, and e is the orbital eccentricity. This yields the 
orbital-period evolution rate 

Pb G i m 

Damour, Gibbons and Taylor showed that the appropriate phenomenological limit on G 
is obtained by 

G SPb 



G 2Pu 
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where 5Pb represents whatever part of the observed orbital period derivative that is not 

From the timing of the binary pulsar PSR 1913+16, a bound is 

(see also 



otherwise explained |]37 
obtained: G/G 



1.0 ±2.3) X IQ-^Vr" i 



38| and p^). However, only for 



the orbits of bodies which have negligible gravitational self-energies, the simplifications 
can be made that Pb/Pb is dominated by —2G/G term. When the effect of the variation 
in the gravitational binding energy induced by a change in G is taken into account, the 



above bound is somewhat weakened depending on the equation of state 40 



3.3 Stars and G 

Since gravity plays an important role in the structure and evolution of a star, a star can be 
a good probe of time variation of G. It can be shown from a simple dimensional analysis 
that the luminosity of a star is proportional to G^ [0. Increasing G is effectively the 



same, via the Poisson equation, as increasing the mass or average density of a star, which 
increases its average mean molecular weight and thus increases the luminosity of a star. 
Since a more luminous star burns more hydrogen, the depth of convection zone is affected 
which is determined directly from observations of solar p-mode (acoustic wave) spectra 



4^ . Helioseismology enables us probe the structure of the solar interior. Comparing 
the p-mode oscillation spectra of varying-G solar models with the solar j9-mode frequency 
observations, a tighter bound on G is obtained: \G/G\ < 1.6 x 10"^^yr~^ . 



The balance between the Fermi degeneracy pressure of a cold electron gas and the 
gravitational force determines the famous Chandrasekhar mass 



l-Ch 



G-'l'm-\ (12) 



where mp is the proton mass. Since Mch sets the mass scale for the late evolutionary 
stage of massive stars, including the formation of neutron stars in core collapse of su- 
pernovae, it is expected that the average neutron mass is given by the Chandrasekhar 
mass. Measurements of neutron star masses and ages over z < 3 ~ 4 yield a bound, 
GjG = (-0.6 ± 2.0) X lO-^V"^ P • 



3.4 Cosmology and G: Big-Bang Nucleosynthesis 

The effect of changing G on the primordial light abundances (especially ^He) is already 
seen in Eq.(^) and Eq.(|^): an increase in G increases the expansion rate of the universe, 
which shifts the freeze-out to an earlier epoch and results in a higher abundance of ^He. 
In terms of the "speed-up factor", ^ = H/Hsbbn, Yp is well fitted by 

Fp ^ 0.244 + 0.074(^2 - 1). (13) 

If Yp was between 0.22 and 0.25, then —0.32 < AG/G < 0.08, which corresponds to 
G/G = (—0.55 ~ 2.2) X 10~^^yr~^. A similar (more conservative) bound was obtained in 
g|: -0.3 < AG/G < 0.4. 
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3.5 Recent Developments on Go 



No laboratory measurements of G/G has been performed recently (see for older labo- 
ratory experiments). This is mainly because the measurements of the present gravitational 
constant Gq itself suffer from systematic uncertainties and have not been performed with 
good precision. 

Recently, Gundlach and Merkowitz measured Gq with a torsion-balance experiment 



in which string-twisting bias was carefully eliminated [47|. The result was a value of 



Gq = (6.674215 ± 0.000092) x lO^^^m^kg'^s"^. Q As the accuracy of the measurements 
improves, it may be possible to place a bound on a present-day variation of G. It is 
important to pursue laboratory measurements of G/G since they are repeatable and 
hence are complementary to astrophysical and geophysical constraints. 



4 A or Dark Energy 

Finally, we briefly comment on the potential variability of the cosmological constant (or 
dark energy) because in the runaway scenario of dilaton or moduli (p, a/a and G/G would 
close to 



4.1 Evidence for A > 

There are two arguments for the presence of dark energy. The first indirect evidence 
comes from the sum rule in cosmology: 

T.^^ = l, (14) 

where Qi = SnGpi/SHQ is the density parameter of the i-th energy component, pi. The 
density parameter of the curvature, Qk, is defined by Qk = —k/a^Hl. Since the current 
obervational data indicate that matter density is much less than the critical density VLm < 
1 and that the Universe is fiat, we are led to conclude that the Universe is dominated by 
dark energy, Q^e = I — — > 0. 

The second evidence for dark energy is from the observational evidence for the accel- 
erating universe f^ : 



I (n^il + + (1 + 3w)nnEil + zf^'+'^A > 0, (15) 



aHl 2 

where w is the equation of state of dark energy, w = Pde/pde- Since distance mea- 
surements to SNIa strongly indicate the Universe is currently accelerating, the Universe 
should be dominated by dark energy with negative pressure {w < 0). We note that an- 
other argument for negative pressure comes from the necessity of the epoch of the matter 
domination. 



^Only recently, the measurement of G with a torsion-strip balance resulted in Gq = (6.67559 ± 
0.00027) X lO^^^m-^kg^^s^^, which is 2 parts in 10^ higher than the result of Gundlach and Merkowitz 
|(48|. Probably the difference is still due to systematic errors hidden in one or both of the measurements. 
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4.2 Supernova and A 

A current bound on the equation of state of dark energy from supernova data isw< —0.6 
|50| . Future observations of high redshift supernovae/galaxies/clusters would pin down 
the bound on w to w < —0.9. The extent of time variation of dark energy density is 
readily seen from the equation of motion: 

^ = -3{l+w)H. (16) 
Pde 

5 Conclusion 

A short account of the experimental constraints on the time variability of the constants 
of nature {a and G)) was given. Since there are some theoretical motivations for the time 
variability of the constants of nature and the implications of it are profound, refining 
these bounds remains important, and continuing searches for the possible time variability 
of the constants of nature should be made. 
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